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Abstract — In this paper we show how to extend the known 
algorithm of nodal analysis in such a way that, in the case 
of circuits without nullors and controlled sources (but allowing 
for both, independent current and voltage sources), the system 
of nodal equations describing the circuit is partitioned into 
one part, where the nodal variables are explicitly given as 
linear combinations of the voltage sources and the voltages of 
certain reference nodes, and another, which contains the node 
variables of these reference nodes only and which moreover 
can be read off directly from the given circuit. Neither do we 
need preparational graph transformations, nor do we need to 
introduce additional current variables (as in MNA). Thus this 
algorithm is more accessible to students, and consequently more 
suitable for classroom presentations. 
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I. Introduction 

It is a well known fact - already taught in most undergraduate 
courses on circuit theory [1] - that, when the node equations 
for a standard Nodal Analysis (NA) of a linear circuit contain- 
ing only admittances and independent current sources are set 
up in matrix form 

where Vn denotes the vector of node voltages v^^ , . . . , 

of those nodes different from some fixed reference node (0), 
the entries of the node -admittance matrix and the node 
current-source vector can be directly read off from the 
circuit itself. I.e., each diagonal term of in position (z, i) 
is given by the sum of admittances incident with the node ©, 
each off-diagonal term in position (z, j), z ^ j, is described by 
the negative of the sum of admittances connecting the nodes 
© and (2); the z-th entry of the vector is the sum of all 
independent currents leaving or entering the node © with a 
plus sign attached only to those currents directed toward the 
node and a minus sign to all the others. 

This article first appeared in: SMACD'04 Proceedings of the International 
Workshop on Symbolic Methods and Applications in Circuit Design, Wro- 
claw, Poland, 23-24 September 2004. Sidney 2004, pp. 113-116. Due to the 
low distribution of these proceedings, the author has decided to make the 
article available to a larger audience through the arXiv. 

The original research which led to this paper was done while the author 
was with the Institut fiir Netzwerktheorie und Schaltungstechnik, Technische 
Universitat Braunschweig, D-38106 Braunschweig, Germany. The preparation 
of this paper and its presentation at the SMACD'04 workshop was not 
supported by any public funding. 

The author's current (as of March 9th, 2009) address is BismarckstraBe 20, 
D-38518 Gifhorn, Germany. Current e-mail: e.gerbracht@web.de 



As is equally well known, while it is easy to extend 
nodal analysis to deal with circuits, which furthermore contain 
voltage controlled current sources or nullors, massive problems 
arise, when voltage sources of any kind have to be taken 
into consideration, as well. Although this obstacle has been 
basically overcome by the invention of the Modified Nodal 
Analysis (MNA), which gives a universal method for any kind 
of linear circuit, with good reasons most teachers of Electrical 
Engineering seem to be very reluctant to confront their stu- 
dents with the MNA-algorithm, especially in undergraduate 
courses. 

Accordingly, several authors (f^,^) have proposed another 
alternative, the so-called Supernodal Analysis (SNA), which 
seems to be more accessible to students and thus has been 
incorporated into existing undergraduate and graduate course^ 
(SHIt!])- Starting with a linear circuit with admittances and 
all kinds of independent sources, the initial set of equations 
can be reduced to a smaller set; these resulting SNA-equations 
again can be described in matrix form as 

YatVat = Jat, 

where vat is a vector of selected node voltages, Yat is a 
matrix of admittances and J at consists of suitably chosen 
linear combinations of the independent sources. Although in 
the literature (ll, Si) instructions are given how to calculate 
the entries of Yat and J at, as far as we know, no algorithm 
has been developed, yet, which in analogy to standard nodal 
analysis allows one to directly read them off from the circuit. 
This paper was written to remedy this situation. 

Remark: Throughout this paper, to keep notation as simple 
as possible, while we freely talk about circuits with admit- 
tances, all the examples will be linear circuits considered in 
the time domain, which besides independent sources consist 
of resistors with positive conductances (symbolized by capital 
letters), controlled sources and/or nullors. The experienced 
reader will know how to generalize the results, which will 
be presented, to other linear circuits containing inductors and 
capacitors. 

The only notational convention we will strictly adhere to is 
using boldface letters for vectors and matrices. 

Without loss of generality (cp. Il], 1.5.3) we demand that 
all circuits under consideration are connected. 

II. Supernodal Analysis of Linear Circuits with 
Independent Sources 

To keep matters simple at the beginning, in this section 
we will only consider circuits without controlled sources and 

^ Since the first publication of this article, the lecture notes | 5] seem not 
to be publicly available anymore. Nevertheless, those contents of the course 
which were relevant for this paper have been incorporated into |6]. 



nullors. 

The basis of our discussion is the concept of a supernode. 
The definition, we will give, slightly differs both from the 
one presented in 0], as well as from the one in chapter 
4.2. This was done to streamline the formulation of the 
"traditional" algorithm and to prepare for our improvements. 

Definition 1: A subcircuit of a given circuit which is con- 
nected, consists only of nodes and (independent) voltage 
sources, and which is maximal with these two propertied is 
called a supernode. 

Let us remark, that by this definition an ordinary node 
which is not incident with any voltage source is regarded as a 
supernode, as well. Furthermore, any supernode defines a cut 
of the circuit. The contractio^ F of a circuit F, obtained by 
contracting all the branches of all supernodes and removal of 
all resulting loops, is a circuit, which by the initial prerequisite 
of this section only contains resistors and independent current 
sources. During the course of this paper we will call F the 
contraction along the supernodes. 

If, moreover, one removes all branches associated to current 
sources from F, the result is the deactivated circuit in the sense 
of Si. 

A. The Algorithm - Proceeding as in Textbooks 

We are now able to adapt the general algorithm of Super- 
nodal Analysis as presented in 12 S Si. and formulate it in 
our terminology. 

Input: a connected circuit F with n + 1 nodes, named 
(O), (I), . . . ,@7 consisting only of independent sources and 
resistors. (* Node (O) will be our global reference node 
(ground/datum), and we set = 0. *) 

Output: 

1) a set of equations for the node voltages 

^® ' * * * ' ^^^^^ completely describes the circuit F. 

2) a subset of node voltages together with a reduced system 
of equations (i.e. equations containing these variables 
only), the solution of which directly leads to the solution 
of the whole system. 

1) Initialization 

a) Assign node voltages. 

b) Identify all supernodes and mark them; let + 1 
be the number of nodes in the contraction along 
the supernodes F. 

c) Within each supernode define one node as the 
local reference node of the supernode (a supernode 
consisting of a single node only is its own reference 
node). The global reference node should be chosen 
to be a local reference nodfl 

2) For each supernode, express the node voltages within in 
terms of the node voltage of its local reference node and 
the values of the voltage sources, it encompasses. 

^I.e. it cannot be enlarged without losing one or the other attribute, 
^in the sense of |9], exercise 11.4.15. 

^Clearly the global reference node need not be chosen before this step of 
the algorithm. 



3) For each supernode, write one KCL equation in terms 
of the voltages of the local reference nodes; leave out, 
as usual, the one for the supernode containing the global 
reference node. 

For those, who are not accustomed, yet, with supernodal 
analysis, we exemplify the workings of the above algorithm: 




The two voltage sources of this circuit give raise to two 
supernodes (a) and (S) consisting of the nodes (I), (2) and 
(0),(4), respectively. While we can freely choose (2) as the 
local reference node for the supernode (a), the local reference 
node of (b) is supposed to be the node @. 

With node voltages being introduced, the internal structure 
of the supernodes implies 




\ 

- ^01 
+ ^02 




(1) 



where we have partially solved this system of equations for 
the voltages of those nodes, which are not reference nodes. 
Kirchhoff equations for supernodes (a) and @ give: 

- + ^3(^0 - v^) + G^{v^ - v^) = 0, 

- + ^2(^(3) - + ^01 = 0. 

Finally, by using ([T]) and collecting all currents and voltages 
resulting from independent sources on the right hand side, we 
are led to the following system of equations for the voltages 
of the local reference nodes: 



Gi + Gs + G4 —Gi 
—Gi Gi + G2 






-^01 



Gs • ^01 
G2 • V02 
(2) 

This reduced system is completely decoupled from ([T]). 
Furthermore its solution together with ([T]) describes all the 
node voltages of our example circuit. 

B. On the results of the SNA algorithm in general 

The preceding example gives rise to some observations, 
which easily generalize to theorems for arbitrary circuits built 
from admittances and independent sources. Let us suppose that 
within each supernode one local reference node has already 
been chosen. 

Theorem 2: If a supernode carries the structure of a tree 
then the voltage of any node within the supernode is uniquely 



expressible as the voltage of the local reference plus the 
sum, relative to the path orientation, of voltages of those 
independent sources along the unique path from the reference 
node to the node under consideration. 

Thus, step[2]of the algorithm can be successfully carried out, 
iff either every supernode is a tree, or any loop of independent 
voltage sources within any supernode satisfies KVL (which 
can be easily checked by identifying a spanning tree within 
each supernode). 

For the sake of completeness, we note that branch voltages 
along admittances which at both ends are incident with the 
same supernode can be determined without any calculations. 

Theorem 3: If one of the above assumptions holds, the 
system of equations resulting from step [3] of the algorithm 
of supernodal analysis can be reduced to 

YtvVat Jat, (3) 

where vat is the vector of voltages of local reference nodes, 
Y AT is the node-admittance matrix of the contraction along the 
supernodes F and J at is a vector of currents induced by the 
independent sources, which will be specified more precisely 
below. 

C. Reading off the Entries from the Circuit 

Since Y^v is the node- admittance matrix of the contraction 
F, there is a simple rule how to fill in the entries just by 
inspection. 

Rule 1: Each element on the main diagonal at position 
(A, A) is the sum of the admittances incident at only one 
enc0 with the supernode (a). The off-diagonal elements at 
position (A, B) is the negative of the sum of those admittances 
connecting supernodes (a) and (S). 

In the future - for obvious reasons - when we talk about 
admittances "incident with a supernode" we will only consider 
those, that are incident at only one end. 

Without loss of generality, let us assume that all supernodes 
of F are trees. By this additional assumption paths within each 
supernode are unique. Thus any admittance Gk connecting 
two supernodes (a) and (b) defines a unique oriented path 
, consisting only of independent voltage sources and the 
admittance Gk itself, as sketched in the figure below: 

(2) . . ® 

The path is uniquely defined as the path which starts in 
the local reference node @ of (S), passes through the node of 
(S) incident with the admittance G^, along Gk to that node of 
(a) incident with the other end of Gk and finally ends in the 
local reference node (a) of (a). When this path is traversed in 
the opposite direction, we will call it 7^^ . Now let be the 
sum, relative to the orientation of this path, of those voltage 

^ As noted above, those admittances incident with only one supernode at 
both ends have to be neglected. 



sources which are part of the path 7^^ , with a voltage being 
counted positive if the voltage source in question and the path 
are oriented in opposite and negative, when the orientations 
are the same. Consequently we have = 0, if G/c directly 
connects two local reference nodes. We are now in a position 
to formulate the fill-in rule for the vector J n - 

Rule 2. Let {im} be the current sources and {Gk} be the 
admittances of the circuit F. Let (a) be a supernode in F, not 
containing the datum node. Then the entry of J at at position 
A is given by 

im-im^(A) ir-ir^(A) G fc : G fc Incident wlth (a) 

(4) 

where im ^ @ means, that the current source Zrn is incident 
with (a) and directed toward (2), and v ^ @ ^^^^ 
incident with (a) but directed away from (a). 

Returning to our first example, in the next figure we see the 
relevant paths 7^2 and 7^3 , which result in the entries G2 • ^02 
and G3 -^01 of the right hand side of equation 




For the convenience of our readers we give a second example 
with the supernodes marked and the relevant "admittance 
paths" sketched. 




Now it is easy to read off the resulting reduced system of 
equations. The interested reader is called upon to check each 
entry by using the original algorithm, to see the advantage of 
our approach. 
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Although this exemplary circuit is highly artificial, it helps 
to stress yet another point, namely that in many circuits any try 
to set up nodal equations with the help of source shifting and 
substitution by Norton equivalent subcircuits is highly difficult 
or nigh impossible. 

D. Circuits with Controlled Sources 

The algorithm of nodal analysis can be easily generalized to 
circuits containing voltage controlled current sources. Clearly 
supernodal analysis and the two rules given above are suitable 
for setting up the equations when circuits with voltage con- 
trolled sources of any kind have to be analyzed. In this case, 
however, the concept of a supernode has to be extended in the 
following way: 

Definition 4: A subcircuit of a circuit which is connected, 
consists only of nodes and independent as well as the con- 
trolled branched of dependent voltage sources, and which is 
maximal with these two properties is called a supernode. 

By temporarily treating controlled voltage sources as inde- 
pendent sources (the "taping" of [4]) and using this modified 
definition, again we can use the above algorithm plus our two 
fill-in rules. To complete the algorithm, in an additional fourth 
step we need to replace the dependent voltages and currents by 
the controlling node voltages (i.e. the untaping of [4]). In this 
more general case, we cannot expect a clean-cut partition of 
the node equations into those, describing the "inner workings" 
of each supernode and those resulting from the contracted 
circuit only. Nevertheless, our rules lead to a simple to set 
up interim result. 

Let us finally consider current controlled sources: When 
supernodes have been introduced, controlling branch currents 
fall into two classes, those inside and those outside a supern- 
ode. Those outside are currents through admittance branches 
and thus can be easily substituted by the corresponding set 
of node voltages as variables. In the second case, with the 
controlling branch being part of a supernode, if again we can 
assume that each supernode is a tree, then by induction we can 
show that at least one end of the branch under consideration 
has to be the root of a rooted tree of voltage sources. The 
leaves of this tree are nodes through which currents induced 
by admittance branches or current sources enter the supernode. 
Now, using KCL at each node of the rooted tree, it is easy 
to describe the controlling branch current as the sum of all of 
these currents. 

^ The reader should keep in mind that a controlled source consists of 
two branches; thus we have to specify which branch should belong to the 
supernode. 



E. Circuits which contain Nullors 

It was the authors initial hope that the approach of supern- 
odal analysis by inspection could somehow be made to work 
on circuits containing nullors as well. This hope was shattered 
by the example belo^^|Z|. 

Though for circuits containing nullors, again there is a 
reduced equation system which is of the form^ArVAr = Jn 
and although there are hints of how to fill out Yat, in general 
there does not seem to be any suitable definition of the notion 
of a supernode which would lead to a short cut, by which 
the tedious setting-up of equations could be avoided: On the 
one hand supernodes in circuits with nullors should contain 
the norators (thus defining Kirchhoff- surfaces for which the 
equations have to be set up); on the other hand the voltages 
of nodes connected by nullators to supernodes (in the sense of 
definition [T]) are known from the beginning as well, and should 
not appear as variables in their own right. As our example 
shows, these two demands do not lead to any well to teach or 
easy to apply fill in rule. 
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-(G1 + G2) G1 + G2 + G5 

— G3 G3 + G4 



-G5 
-G4 



-G2 ■ vol ^ {Gi ^ G2) ■ V02\ 
-G3 • ^01 + (G3 + G4) • V02 J 

Fortunately, in case of admittance circuits with nullors, 
supermesh analysis or RLA, as it was presented in algorithm 
3.2 and chapter 5 of can be shown to be universal, and 
thus can be used for all circuits without cutsets of independent 
current sources and nullators. 

III. Conclusion 

Past experience in Braunschweig with an undergraduate 
course of about 150 students |7] has shown that the modi- 
fied SNA-algorithm, with fill-in rules, as presented above, is 
quite easy to teach. Moreover, the students naturally adapted 

^The reader should note the particular choice of nodes (T) and (4) as local 
reference nodes, which from a practical point of view clearly is absurd, but 
from an algorithmic point of view is a distinct possibility. Furthermore we 
have refrained from setting = 0. 



to it and after only a small amount of training applied it 
successfully to a number of problems. Thus it seems that this 
algorithm together with its counterpart, the supermesh analysis 
(for nullor circuits), has reached a state of maturity, that makes 
it perfectly suitable for presentation as "universal tools" in any 
course on circuit theory. 
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Note added to the Electronic Version 

In this electronic document, some small typographical errors 
of the printed version were corrected. The figures, though still 
drawn by hand, now are from the original coloured drawings. 

Furthermore, for the convenience of the reader the ab- 
stract has been rewritten, and keywords, ACM and MSG 
classifications, and a short CV according to IEEE stan- 
dards have been added. One URL has been added to the 
bibUography. (March 9th,2009) 



